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ABSTRACT 

In this paper the expected properties of the velocity fields of strongly lensed arcs 
behind galaxy clusters are investigated. The velocity profile along typical lensed arcs 
is determined by ray tracing light rays from a model source galaxy through parametric 
cluster toy-models consisting of individual galaxies embedded in a dark cluster halo. 
We find that strongly lensed arcs of high redshift galaxies show complex velocity 
structures that are sensitive to the details of the mass distribution within the cluster, 
in particular at small scales. From fits to the simulated imaging and kinematic data 
we demonstrate that reconstruction of the source velocity field is in principle feasible: 
two dimensional kinematic information obtained with Integral Field Units (IFU's) 
on large ground based telescopes in combination with adaptive optics will allow the 
reconstruction of rotation curves of lensed high redshift galaxies. This makes it possible 
to determine the mass-to-light ratios of galaxies at redshifts z > 1 out to about 2-3 
scale lengths with better than ~ 20% accuracy. We also discuss the possibilities of using 
two dimensional kinematic information along the arcs to give additional constraints 
on the cluster lens mass models. 

Key words: gravitational lensing - techniques: interferometric - galaxies: high red- 
shift - galaxies: kinematics and dynamics - cosmology: dark matter 
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1 INTRODUCTION 

In recent years, the high redshift universe has become a fo- 
cus of attention for observational cosmology as a testbed 
for theoretical models of galaxy formation and cosmology. 
The detection of high redshift galaxies and quasars has 
led to the first observat i onal tests of reio nisation models 
(lHaiman fc Holder! 1200 A ICiardi fc Ferraral feOOffl and high 
redshift absorption line systems have helped to constrain 
spin temperatures an d the fraction of cold neutra l gas in 
high redshift galaxies (iKanekar fc Chengaluil l2003l bource 
counts in the sub- mm and the infrar ed have constrained 
the star formatio n at high redshifts (jHughes et al J 11998: 
Blain et al. 1999). However, despite these tremendous ad- 
vances in the field, a few questions pertaining to the high 
redshift universe and the evolution of galaxies remain diffi- 
cult to address. Foremost, there is the open question as to 
the role of the dark matter in galaxy evolution. Even in the 
local universe, the presence of dark matter can only be in- 
ferred indirectly through its gravitational effect; for example 
by measuring the rotation curves of galaxies. 

Rotation curves have now been measured for a 
large sample of nearby galaxies, both in optical wave- 
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lengths (iMathewson et al.1 Il992t IPersic fc Saluccil \E 
Palu nas fc Williamsl l200(t) and using the 21cm emission 
line o f neutral hydrogen (Verheiien & Sancis 3 l2Q0lHSwater3 
1999). They have improved our knowledge of the systematics 
of dark matter in nearby galaxies, but they have also led to 
a number of new questions that still need to be addressed in 
theoretical models of galaxy formation. The strong depen- 
dence of rotation curve shape and amplitude on the optical 
luminosity indicates a tight coupling between luminous and 
dark matter that is not expected from the current models 
(IPersic et al.lll996l : iDonato et alJfeOOA . There are also in- 
dications from the rotation curves of low surface brightness 
galaxies that galactic dark matter haloes contain a constant 
density c ore, as opposed to the r -1 cusp predicted b y sim- 
ulations (Ide Blok et al.ll200ll : Ide Blok fc Bosmall2002F) . 

Using Ha emission, rotation curves of galaxies have 
been measured accurately up to redshifts of about 0.3. Ha 
emission has been detected out to much larger distances 
and has been used to constrain th e star formation rat es in 
galaxies at redshifts greater than 2 (iBunker et al.lll999h . At- 
tempts have also been made to meas ure the rotation curves 
of galaxies up to redshifts of z ~ 1 (IVogt et alJll996l . [l£ 



Hu dson et alll998tlB5hm et all20oA . but at these redshifts 
it is difficult to determine the inclination angle of the rotat- 
ing disc accurately and to obtain sufficient spatial resolu- 
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tion. So far the limited spatial resolution has only allowed 
determinations of the total velocity dispersion rather than 
actual rotation curves: therefore most authors concentrated 
on constraining the evolution of the Tully-Fisher relation 
iTiillvfeFiaheJlflri). 

Galaxies that are at a higher redshift but are strongly 
lensed may well be magnifie d by factors of 10 or more (e.g. 
Bla ndford Sz Naravanl Il99^ ). Exploiting this effect, gravi- 
tational lensing has already been successfully used as a 
tool to study high redshift s ources in greater detail than 
would be otherwise poss ible (jPello et alJ Il999l : lEllis et all 
l200ll : lRichard et all200fl . Lensed galaxies may have a large 
enough area and luminosity to make a measurement of 
their rotation curve possib le, despite their larger distance. 
iNarasimha fc Chitrd (|1993i) demonstrated that it is in prin- 
ciple possible to measure the kinematics of galaxies that are 
much more distant in this way. Since integral field spec- 
trometers had not been developed at the time of the publi- 
cation of that paper, it focused on measuring the velocities 
along straight arcs. Long slit spec troscopy of straight arcs 
has been carried out successfully by Bun ker et alJ $2000) and 
Meh lert et alJ ([200 ill . Recently, improvements in resolution 
and sensitivity of integral field spectrometers have made it 
possible to obtai n kinematics of mode rately magnified back- 
ground galaxies (|Swinbank et ajj 2003). 

At intermediate redshifts between z ~ 0.3 — 1, gravita- 
tional lensing also provides constraints on the dark mat- 
ter content of the inner 1 — 200 kpc of lens ga laxies 
(e.g. iKoopmans et alJ Il99fifc E ooomans and 

lens clusters (e.g . ISmail Dickinson! Il995l : iMellierl [l 

iKneib et alJ 1200*3) . Determinations of the dark matter dis- 
tribution using gravitational lensing rely on accurate data to 
provide sufficient constraints on the lens mass model. Degen- 
eracies in the lens model usually exist and may have impor- 
tant consequences for determinations of cosmol ogical param- 



eters and mass density profiles from lensing IWilliams e^al 
)inl l2003l:lMeneghetti et al.ll2004l : i Dalal et al 



). Additional constraints help to break such degenera- 
cies. For galaxy lenses, spectroscopic studies have been used 
to break degeneracies in t he lens models and constr ain their 
mass distribution further (|Koopmans &; Treul l2002'). 

In this paper we investigate the possibilities of measur- 
ing rotation curves of lensed galaxies with current and up- 
coming instruments and discuss the additional constraints 
on the foreground lens mass model that may be obtained 
from two dimensional spectroscopic data. 

We begin by describing the theoretical framework, the 
ray-tracing method and the source and cluster lens models in 
§[3 In §0 we briefly discuss the optical properties of known 
lensed arcs. In §|I| we present the kinematics of simulated 
arcs and discuss their generic properties. In §|5| we discuss 
how the kinematic profile of the source can be reconstructed 
from lens modeling and discuss their dependence on the lens 
mass model. We address the observational possibilities using 
current and future instruments - in particular the use of 
Integral Field Units (IFUs) - in §|S| We conclude with a 
discussion and summary in §0 

Throughout this paper we use a standard ACDM 
cosmology with Q\ — 0.7, Q m — 0.3 and Ho — 
70 krns" 1 Mpc -1 . 



2 SIMULATING DATA CUBES OF LENSED 
HIGH REDSHIFT GALAXIES 

2.1 Lensing theory 

A background source at a redshift z s that is located at a 
position (3 will appear at a position on the sky if it is 
lensed by a massive foreground object at redshift z\ so that 



(1) 



where Dls and Dos are the angular diameter distances from 
lens to source and from observer to source, respectively. The 
deflection angle a at position on the lens plane is given by 
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where k is a dimensionless quantity related to the surface 
mass density £ through, 



K (0) = £(0) x 
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The magnification of an image at position 6 of a point 
source at position (3 is given by 
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[1 - k(6)] 2 -7 2 (6>)' 



(4) 



where 7 = (j x ,J y ) is the total shear at the image position, 
which is related to the deflection angle, and hence the lensing 
mass distribution, through: 
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For extended sources, different parts of the source will be 
magnified by different amounts, leading to differential mag- 
nification (|Blainl ll99£0. Given a magnification map on the 
source plane /jl s (/3) and the source flux F s (/3), the observed 
total flux is 



F a (p)^ s ((3)d 2 p. 



(7) 



The type of observations we are interested in here, 
namely spatial spectroscopy using integral field units, are 
characterised by several channels at frequencies z/, of band- 
width Aza The observed flux in an interval between v and 
v + dv would be 



(8) 



F{y)dv= J F a (v,p)fjL a (f3) d 2 f3dis, 
and so, for finite band widths, 



F(v) 



-n 



F a (u,/3)ns (13) d'Pdv. 



(9) 



Thus, if /i s (/3) is known from the lens model, it is straight- 
forward to obtain the total observed flux at a given wave- 
length, given a model of the source. This can be done simply 
by summing over all the pixels of a map that is the product 
of the source flux at the observed wavelength and the source 
magnification map. 
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2.2 Modeling the source and creating mock 
kinematic data 

In order to create a model data cube for the source, we 
need to make assumptions about the gas distribution and 
its kinematics. 

For the l atter, we take the 'Universal Rotation Curve' 
as derived bv lPersic et al.l (|l99(t ). We assume an L* galaxy, 
for which their equations reduce to: 



1.4a; 1 - 2 



2 + 0.6) 1 - 4 



+ 



l.lx 2 
z 2 + 2.2' 



(10) 



with x = R/Ropt = R/3.2h being the radius expressed in 
units of i?opt, which is the radius encompassing 83% of the 
total light; h is the disc scale length. We set the parameter 
Vq = 201.5 km/s. Defining the radius expressed in disc scale 
lengths, Xh = R/h, we get: 



V*(x h ) = V ' 



O.lxl 



(O.lxl + 0.6) 1 - 4 O.lxl +2.2' 



(11) 



In reality, galaxies do not follow the Universal Rota- 
tion Curve exactly, and galaxies with different masses will 
have differen t shap es for their rotation curves. Furthermore, 
IPersic et al.l (|l99ft) derived their equations from rotation 
curves of galaxies in the local universe, and their results 
may not be directly applicable to the high-redshift galaxies 
we are studying here. However, the exact shape and ampli- 
tude of the rotation curve is not critical for the study we 
present here and different assumptions would not affect our 
conclusions. 

For the gas distribution, we assume an exponential ra- 
dial profile: 



I = I exp(-Xh) 



(12) 



We adopt a radial scale-length of h = 0.4arcsec, cor- 
responding to approximately 4 kpc at the assumed source 
redshift of z s = 1.5. Furthermore we assume that the disc of 
the galaxy is inclined at an angle of 50 degrees with respect 
to the line of sight. The vertical density distribution of the 
gas is assumed to be exponential as well, with scale-height 
zo equal to l/20th of the radial scale- length h. This means 
we effectively assume the galaxy to be razor thin, but the 
effect of a larger vertical scale-height on the simulated data 
cube is very small for the moderate inclination of 50 degrees 
we assume here. 

To create the model data cube, we use the GALMOD 
task in GIPSY: The Groningen Image Processing SYstem. It 
uses a Monte- Carlo integration to fill the model galaxy with 
small gas clouds, following the exponential distribution given 
above. For each cloud, the radial velocity is calculated on 
the basis of its position in the disc and the rotation velocity 
at its radius, and each cloud is assumed to emit a gaussian 
emission line with an intrinsic velocity width of 10 km/s. The 
final data cube consists of channel-maps spaced ~ 8 km/s, 
with spatial pixels of 0.02arcsec, or h/20. 

In Fig.[U we show a cut through the data cube along 
the major axis of the galaxy. 

From the data cube, we derive an image of the in- 
tegrated gas emission by adding up, at each pixel posi- 
tion, the signal in the individual channel maps. The veloc- 
ity field is derived by fitting gaussians to the line profiles. 
It is shown in the inset in Fig.Q and resembles those ob- 
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distance along slit (arcsec) 

Figure 1. Position- velocity diagram of the source along the major 
axis. Contours correspond to the maximum intensities at 1, 3 and 
5 radial scale lengths. The green line indicates the rotation curve, 
multiplied by sin(50°). The inset in the top left shows the velocity 
map of the source. The distances along the axes are in arcseconds. 
Different colours indicate radial velocities, ranging from -175 to 
175 km/s, as indicated by the colour bar; the contours range from 
-150 to 150 km/s with intervals of 25 km/s. 



served in nearby sp iral galaxies (|Yerheijen |5ancisj[2QQll 
iGarrido et al.ll200l. 

In the first instance, we use a very high number of small 
gas clouds in the Monte- Carlo integration. This results in a 
very smooth gas distribution and velocity field. We can sim- 
ulate a very patchy gas distribution by limiting the number 
of clouds to a very small number. We will discuss how such 
strong substructure in the source will affect our results in 
5 4.2 and 8 5.5. 



2.3 Modeling the lensing potential 

Simulations of strongly lensed ar cs in clusters have 
been p erformed by severa l groups (|Bartejmaim fcWej^ 
I l994l: IWu fc Maol Il99d : Irlamana fc Futamasel Il997|: 
Bartelmanri^^aL| | l99d jW ambsga nss et alJ [1998; 
iMeneghetti et all l200ll . 12003) . Most of these studies 
make use of a cluster model from N-body simulations. In 
contrast, the basic approach in our work is to use an ana- 
lytic, parametric cluster model. Analytic models have been 
used to predict th e arc stati s tics fo r diff erent cosmologies b y 
IWu fc Maol (ll99(t ). ICooral dl999l) and iQguri et alJ fcOOlft . 
finding a much lower cross section for the formation of 
arcs than when clusters are modele d directly from N -body 
simulations. As pointed out by IMeneghetti et all (|2003) 
massive cluster substructure and assym metry are important 
and e xplain part of the discrepancy. IWambsganss et aTl 
([2004T) recently demonstrated that the source redshift is 
also an important factor affecting the statistical incidence 
of strongly lensed arcs. Furthermore, since small scale struc- 
ture may have a very significant, local effect, the actual 
structure and shape of the arcs will be influenced strongly 
by the presence of nearby mass concentrations. A s shown 
bv IMeneghetti et al.1 (|200(tl and lFlores et al.1 fcOOGb cluster 
galaxies are unlikely to be massive enough to affect the arc 
statistics - more massive substructures and asymmetries 
are needed. However, structures on galaxy scales are very 
important for detailed modelin g of observed arc s ystem s 
(iKassiola. Kovner Ford llOol iBroadhurst et~aH Hf" 
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When such structures lie close to extended arcs, the shape 
of the arcs may be strongly affected. Therefore, a realistic 
parametric cluster model that includes the individual 
galaxies is used here, incorporating the substructure that 
can crucially affect the appearance of strongly lensed arcs 
behind cluster lenses. 

In order to simulate the lensing potential of a clus- 
ter, we create 6 realizations of a mock cluster at redshift 
z\ — 0.3 by arranging 70 galaxies randomly around the 
centre of a common dark matter halo. The distribution of 
galaxy positions follows a gaussian with a standard devi- 



ation of var» 



65 arcsec around the cluster centre. The 



angular distribution is isotropic. Most clusters are expected 
to have a moderately elliptical potential and galaxy distri- 
bution. Such an ellipticity introduces additional parameters 
into our model but does not significantly affect our main 
results. We will discuss the effect of lens ellipticity and pro- 
file in detail in $ 15.41 Each galaxy is modeled as a singular 
isothermal sphere (SIS) of the form, 



E(r) 



2Gr' 



(13) 



For each galaxy, the velocity dispersion cr v is determined ran- 
domly from a gaussian distribution of mean a® = 130 km s _1 
and standard deviation var° = 92kms _1 . In order to inves- 
tigate the effect of the amount of substructure due to the 
cluster galaxies, we parameterise the relative contribution 
of the cluster galaxies to the total cluster mass, using a sin- 
gle parameter, T. The mean velocity dispersion for the SIS 
galaxies is a v — Ht^, where a° — 130 km s -1 . 

In addition to the galaxies, we add a cored isothermal 
sphere (CIS) at the centre of the cluster, to model both the 
effect of a common halo and a central cD galaxy. The surface 
mass density of the CIS is given by 



E(r) 



2G 



(14) 



where r c is the core radius and <j c \ is the velocity dispersion 
at r — > oo. The CIS is centred on the mean galaxy position. 
Its maximal velocity dispersion is set to cr™ ax = 1281 kins" 1 , 
and it has a core radius of 10 arcsec, corresponding to 
110 kpc at z = 0.3. For values of F > 0, that is, for an 
increased mass fraction in galaxies, the cluster velocity dis- 
persion is decreased below cr™ ax to keep the total mass in- 
side the Einstein radius constant at Mer = 1.7 x 10 14 M©. 
By varying the relative values of T and the cluster veloc- 
ity dispersion a c \, the degree of 'dumpiness' in the clus- 
ter is changed while keeping the Einstein radius constant 
at 32 arcsec. We set the size of the simulated strong lens- 
ing image to 2' x 2', which corresponds roughly to the field 
of view of the Hubble Space Telescope (HST). The surface 
mass density maps of 6 different model clusters are shown in 
Fig. El together with the simulated arcs (contours). The dif- 
ferent cluster models are summarised in Table 1. Note that, 
for all models, the individual galaxy positions and redshifts 
remain fixed. 

Our particular choice for the structure of the cluster 
model is simple. It is motivated mainly by the requirement 
that the total mass distribution of the clus ter is close to 
the N avarro, Frenk and White profile fNFW. lNavarro et alJ 
1997). The galaxy and cluster halo mass distribution used 



Model 


r 






var CT 


A 


1.0 





131.3 


91.7 


B 


0.9 


558.5 


118.2 


82.6 


C 


0.7 


914.9 


91.9 


64.2 


D 


0.5 


1109.5 


65.6 


45.9 


E 


0.2 


1255.3 


26.2 


18.3 


F 


0.0 


1281.2 


0.0 


0.0 



Table 1. Parameters for 6 simulated mass models of clusters. The 
relative contribution of the individual galaxies, parameterised by 
T is listed together with the central velocity dispersion of the 
smooth cluster halo. The Einstein radius is 32 arcsec for all mod- 
els. The tabulated values for a c \, cr v and var^ are in kms -1 . 



here gives a total mass profile that is nearly identical to an 
NFW profile. 

In order to quantify the general lensing properties of 
the lensing cluster models we calculate magnification and 
surface mass density maps. These maps illustrate important 
general lensing properties of the cluster lens models. They 
are calculated using gLens by mapping small triangles from 
image plane to source plane as described in Moller & Blain 
(1998). Figs.|2 and E| show the surface mass density of the 
clusters we used as mass model template and the resulting 
magnification patterns on the image plane. The magnifica- 
tion maps show the magnification \i of a point source as a 
function of position. The corresponding magnification maps 
on the source plane are shown in Fig.0] The high magnifica- 
tion regions on the image and source planes are the 'critical' 
and 'caustic' lines respectively. Note that the structure of 
the critical lines, for Y > 0.5 are very similar to the ones ob- 
tained from lens mass reconstructions of observed clusters 
jKneib et al JllflflS iBroadhurst et aTlEooi . 

These figures demonstrate that increasing the mass in 
the individual galaxies, and thereby increasing the amount 
of substructure, creates more strongly distorted critical lines. 
Any long arcs of highly magnified sources that are produced 
along the strongly curved sections of the critical lines will 
appear broken and distorted. Therefore, the probability of 
observing broken and distorted arcs increases with increas- 
ing fraction of the total cluster mass that resides in indi- 
vidual galaxies. For most of the remainder of this paper we 
will use a model with a very strong amount of substruc- 
ture (model A) . This model probably represents an extreme 
case. Since we expect the accuracy with which the source 
can be reconstructed to decrease with increasing amount of 
substructure, this will be a 'worst-case' scenario for any re- 
construction method. 



2.4 Generating a lensed data cube using 
ray-tracing 

In order to obtain simulated images of the source at different 
wavelengths, we use the ray-tra cing code gLens. previously 
used for several le nsing studies (Mol ler Sz Blainlll998l .[2001: 
Moll er et al1 l2003). All pixels of a given image are mapped 
from the image plane onto the source plane. The flux of 
the zth pixel in the image plane, fi = F I (6 
flux on the source plane F s from F l {6i) 



, is set using the 

S (a. ^.£>ls\ 



but is not the 



This method works well and is very robust 
most efficient way to calculate images of extended sources 
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Figure 2. The surface mass density and lensed arcs of the model 
lensing clusters. The panels are for models A-F in Table 1. The 
colour-scale shows the surface mass density of the cluster. The 
lines show the contours that contain 10, 68 and 90 per cent of 
the total source flux. The value of T decreases from left to right 
and top to bottom. Note that for larger values of T, that is, for 
larger amounts of substructure, the arcs are broken and distorted 
in several places due to the presence of the individual galaxies. 



For example, many pixels are mapped to the empty re- 
gions in the source plane. Ray-tracing of these 'empty' rays 
could be avoided usi ng some sort of adaptive algorithm 
(Mol ler fc Blainl l200lT) . However, for the purpose here the 
non- adaptive approach is sufficient. The ray-tracing tech- 
nique is extremely accurate. Numerical artefacts appear only 
when there is a strong mismatch between the source and 
image plane resolutions or dimensions. In this paper, the 
image plane has a size of 120 x 120 arcsec and a resolution 
of N = 1200 x 1200 pixels. This image plane is mapped 
onto a source plane of 800 x 800 pixel resolution, covering 
an area of 40 x 40". With these settings numerical errors 
are negligible. 



3 PROPERTIES OF LENSED ARCS 

Seve ral strongly lensed arcs have bee n discovered to date 
fe.g.lFort et al|l988l:lKneib et al.ll995l: 
bladders et al.l 120031 : iBroadhurst et al. 



Luppino et all 19991: 



2005). The arcs in 
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Abell 2218 and Abell 370 are perhaps the most notable 
of these, being up to 20 arcsec long and 2-3 arcsec wide. 



Figure 3. The image magnification maps of the model lensing 
clusters. The panels are for models A-F in Table 1 and show the 
magnification of point sources as a function of image position. 



Smaller arcs have been discovered in some other clusters, 
like CL 1358+62. In these clusters, the arcs are usually less 
extended in both directions with lengths of a few arcseconds 
and widths ~ 1 arcsecond. It is important to make a distinc- 
tion here between arcs that are produced by galaxy lenses 
and arcs that are lensed by massive clusters. The lengths and 
widths are quite different. An example of an arc produced 
in strong galaxy-galaxy l ensing is the most p rominent arc 
in the Ultra Deep Field (iBlakeslee et al.ll2004l UDF). The 
lensing galaxy is a field elliptical, much less massive than 
the central cD galaxies in massive clusters. The arc is only 
~ 0.3 arcsec wide and about ~ 1 arcsec long. 

We show the resulting images of our simulated lensed 
arcs as contours in Fig.|5| The simulated arcs for models A 
and B are very similar to those actually observed for several 
cluster lenses and show the main features clearly: a broken 
structure caused by the individual galaxies in the lens. The 
arc lengths and widths of ~ 10 arcsec and ~ 2 arcsec, respec- 
tively, are very comparable to observed widths and lengths. 
The smoother potentials of models C-F lead to continuous, 
smoother arcs. The appearance of the lensed arcs is very 
similar to observed arc systems for models A and B and we 
therefore use only these two models for the remainder of this 
paper. 
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Figure 4. The magnification map on the source plane for the 
model lensing clusters. As for Fig-EI the panels are for models 
A-F in Table 1. The scale is logarithmic for all plots. 

4 THE VELOCITY FIELDS OF SIMULATED 
ARCS 

4.1 Smooth sources 

Using the ray-tracing procedure described in $ 12.41 we cal- 
culated the individual channel maps of the lensed source, 
modeled as described in $ 12.31 The velocity fields of the 
strongly lensed arcs are determined by fitting gaussians to 
the line profiles. In this section, we do not include obser- 
vational effects, like seeing, coarse instrument resolution or 
noise. These will be discussed in §|H] 

We show the lensed velocity fields of the main arcs for 
models A and B in Fig.|5] The resulting velocity structure 
is very complex. Gravitational lensing into multiple images 
produces a very distorted and asymmetric velocity structure 
along the arc. When velocity information of a lensed arc is 
available, the source can be essentially broken down into sev- 
eral smaller components which cover different parts of the 
source plane and are therefore all magnified and distorted 
in a different way. This effect of differential magnification 
was a l ready discussed in $ 12. H and, in a different context, by 
Blain] (jl999h . The small difference between models A and 
B in terms of the mass distribution within the cluster (in 
model A, the galaxy that is closest to the arc has a velocity 
dispersion of 248 km sec -1 , in model B it is 223 km sec -1 ) 
translates into a noticeable and measurable difference in the 
velocity fields: in model B, the regions with approaching ve- 
locities are more strongly magnified relative to the receding 
side. The left and right panels in Fig.|H|show the position ve- 
locity diagrams for model A and model B, respectively. The 



positions of the curved slits used to generate the position 
velocity diagrams are indicated in Fig.|5]by the white lines. 
They have a width of 1 pixel or 0.3 arcsec. A broader slit 
would increase the velocity range at each position element 
and the curves would become broader. The two lower sets 
of panels show the results if the original slit is displaced by 
±3 pixels or ±1 arcsec. This changes the regions along the 
velocity field of the lensed arc that are probed, leading to 
changes in the shape of the curve. In particular, remarkable 
features like strong asymmetry in the position-velocity dia- 
gram may result; as seen in panels (c) and (d). Note that 
none of the lensed position-velocity slices resemble those 
of unlensed galaxies. It is also noteworthy that the lensed 
position- vel ocity slices shown h ere look very similar to those 
observed bv lPello et all (Il99ll Fig.4). 

There is another simple and drastic effect that grav- 
itational lensing has on the velocity structure of sources, 
which can be demonstrated easily by looking at the total 
flux emitted as a function of channel. Fig. shows the flux 
ratio F(y)/Ftot in the main arc (top) and the correspond- 
ing ratio F(y)/Ftot for the counter arc (bottom). In both 
cases, F(y) and F to t are given by eqs. 9 and 7, respectively. 
Comparison with the unlensed case, as shown by the dotted 
line in Fig.[7| shows the effect of differential magnification 
very clearly. For most models, the parts of the source with 
velocities around — 100 km s -1 are magnified more strongly 
than the rest of the galaxy. For a symmetric galaxy such 
a profile is a strong indication of differential magnification. 
For a single arc, this effect may also be reproduced without 
differential magnification w hen there is strong asym metrical 
substructure in the source (Richter & Sancisi 1994). Differ- 
ential magnification can either enhance or partially cancel 
intrinsic asymmetries in the source. However, since the dif- 
ferential magnification is strongly influenced by the small 
scale - and hence local - structure in the lensing poten- 
tial, lensing will in general produce different asymmetries in 
different arcs of the same source. In the given case, a com- 
parison between the fluxes for the main arc, in the top panel 
of Fig.[7| and the counter arc, in the bottom panel, shows 
that the asymmetries are induced by lensing. In this way, 
asymmetries induced by lensing can be distinguished from 
intrinsic asymmetries, which affect all arcs of a given source. 

4.2 Source substructure 

Throughout we have assumed a smooth light distribution 
for our source galaxy. At high redshifts the merger rate 
is expected to be m uch higher than in the local universe 
JPatton et alJ f2002V and consequently galaxies may show 
far mo re substructure a nd kinematic signatures of merger 
events (iNaab fc Bur kert 2003). Currently, there is still some 
uncertainty as to how much of the difference in the morphol- 
ogy and observed light distribution of high redshift galaxies 
is due to intrinsic differences and how much is due to the 
fact that the optical at a redshift of z ~ 1 corresponds to 
restframe UV. Due to regions of star formation, local late- 
type galaxies are found to show much more substructure 
when observed in the UV. How would our results change 
if the source light distribution is intrinsically less smooth ? 
As long as the underlying kinematics does not change, we 
find that the dumpiness of the source has little effect on the 
results discussed in the previous sections. This is illustrated 
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Figure 5. Velocity fields for an extended source lensed by a galaxy cluster. The parameters for the lens and source are as in the previous 
figure for model A and B. The white lines indicate slices along which the posit ion- velocity diagrams are determined (shown in Fig. 6). 
The colour scale is the same as for the inset in Fig.Q] 
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Figure 6. Position-velocity diagram along the main lensed arcs. Only the velocity profiles for models A and B are shown. Panels (a) 
and (b) show the position- velocity diagram along a line following the brightest pixels on the arc (the thick white line in Fig.|5|. In panels 
(c) and (d) this line is offset by +3 pixels (« 1 arcsec) in x, in panels (e) and (f) it is offset by -3 pixels in x. The contours correspond 
to the maximum intensities at 1, 3 and 5 radial scalelengths in the input model (cf. Fig.Q. For all offsets the position- velocity diagram 
is heavily distorted with respect to the unlensed case. 
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Main arc system 
(arcs in top right 
of fig. 5) 

Increasing line thickness: 

r=o.5, r=o.7, r=i.o 

Dotted: unlensed 




Counter arc 
(arc in bottom left 
of fig. 5) 
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Figure 7. The flux of the simulated image systems as a function 
of the channel. The ratio of the flux in each channel to the total 
flux in the arc is shown. The top panel shows the results for 
the main lensing arc and the bottom panel shows the results for 
the counter arc. With increasing line thickness, the curves are for 
T =0.5, 0.7 and 1. The dotted curves show the unlensed case. Note 
that, to increase clarity, the dotted curve is normalised differently 
than the solid curves. 



in Fig.|H| where we show the clumpy source and the position 
velocity diagram along the lensed arc, which has the same 
shape and position as for the smooth source model. Apart 
from many discontinuities in the diagram the shape is the 
same as in Fig.|H] However, we still assume dynamically sta- 
ble rotation everywhere within the source. In particular our 
model assumes that there are no major in- or outflows re- 
lated to the source. If these are present, as may be expected 
for a fraction of sources at high redshifts, these will show 
up as clear signatures in the reconstructed velocity profiles 
of the lensed arcs. We discuss in more detail how well sub- 
structured sources can be reconstructed in 5 5.5. 



5 RECONSTRUCTING THE SOURCE 
KINEMATICS AND THE LENS MASS 
DISTRIBUTION 

5.1 Reconstruction algorithm 

In the previous sections we discussed the velocity structure 
of lensed arcs of extended sources. Using a simple model 
for the kinematic structure of the source, we have shown 




distance along slit (arcsec) 

Figure 8. Position- velocity diagram for a lensed clumpy source. 
The source and lens model parameters are as for Fig.|(J model A, 
except that the source now has a clumpy structure as indicated 
in the inset. 



that the observed velocity structure in the arc is complex 
and depends on the mass distribution in the strong lensing 
cluster - especially the mass in the galaxies close to the arc. 

A remaining question is how well the original veloc- 
ity structure of the source can be reconstructed from the 
available channel information of the lensed arc. Clusters of 
gala xies have been mod eled from strong lensing in the past 
(e.g. iKneib et al1 ll998V deriving the cluster potential from 
the positions and shapes of the arcs. This approach works 
very well when the overall cluster potential is to be deter- 
mined. When the source itself is to be reconstructed, it is of- 
ten necessar y to perform more elaborate fits including pixel 
information (iTvson et al.lll998h . 

Using our simulated images for cluster model A, we at- 
tempt to reconstruct the cluster and source model param- 
eters from the image data alone. We assume that galaxy 
redshifts and positions in the cluster are known, but not 
their masses. A number of 14 galaxies are included in the 
cluster model, in addition to a halo of unknown position 
and mass. The free source parameters are position, expo- 
nential scale length, total flux, position angle and axis ratio. 
In total we therefore have 14 + 3 + 6 = 23 free parameters. 
We use a total of Af p i x = 500 pixels selected randomly from 
the region around the arc on the image as constraints. Each 
pixel is mapped onto the source plane and the flux at the 
source plane position is compared with the flux of the source 
model. The positions of the pixels themselves also provide 
a constraint: bright pixels should be clustered more than 
faint pixels. We include this constraint by first calculating 
the flux weighed centre of the mapped source pixels, 



x 



E^pix f 
i=l Jt 



(15) 



and then calculating the distance of each pixel with respect 
to this flux weighed centre. The total 'goodness', £ is then 
calculated as: 



i=l 



xo| 2 /f + {fi 



(16) 



In this equation, vra is the model flux at pixel i. Note that 
we include a dependence on the pixel positions fi in the first 
term, since in our source model brighter pixels are required 
to be more compact than fainter pixels. Also, note that this 
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Figure 9. The colour scales in the top 3 rows show the velocity fields of the 9 best-fitting reconstructions of the smooth source, 
corresponding to lens model A. The bottom row shows the velocity field of the three best reconstructions for the substructured source 
described in $ 14.21 The lines show the isocontours of the corresponding images at 1, 3 and 5 scale lengths. 



definition is only useful for determining the best fit models 
- a meaningful % 2 value can only be defined on the image 
plane. This is done below in § 5.2. 

In order to obtain an acceptable fit, and also to in- 
clude possible degeneracies, we perform the fitting using a 
modified simulated annealing technique, with slow cooling. 
A population of 400 model clusters, initially randomly sam- 
pling the parameter space in a uniform manner, is slowly ad- 
justed in 4000 steps. At each step, a new point in parameter 
space is chosen, sampling the logarithmic parameter space 
using a Monte Carlo Markov Chain (MCMC) method. After 
each step, the new set of model parameters replaces the old 
one with a probability given by 

Preplace = mill ( e <*old-X„ew)/T ^ _ 

where T is the current 'temperature' of the system, which 
is cooled from T sta rt = 10 to T end = 0.001 in the 4000 steps 
logarithmically. From the final sample of 400 models, we 
select the 9 best-fitting models. 



5.2 The reconstructed sources 

The contours in Fig.|5]show isophotes of the resulting source 
reconstructions for the 9 best-fitting lens models. These 
source reconstructions are calculated by ray-tracing the orig- 
inal, 'observed' image to the source plane through the cor- 
responding reconstructed model of the lens. All sources are 



compact and in rough agreement with the input position 
angle and inclination. Most importantly all the best-fitting 
models give source models with an integrated light profile 
that is consistent with an exponential disc profile. 
We define a goodness of fit in the image plane as 

N 2 
i=l 

where 

fi + rrii ( . 

w * = f , > ( 19 ) 

and a — 5 x 10 -5 is our assumed surface brightness error 
in units of the total flux of the source. With this defini- 
tion, x 2 ~ 1 — 3 in all of these cases. The velocity fields of 
the best-fitting reconstructed sources are displayed with the 
colour scales in Fig.|^] Most of the reconstructed velocity 
fields show the same global shape as the input model, but in 
several cases distortions are present, especially in the outer 
parts. These distortions are inconsistent with dynamically 
stable rotation and can be used to distinguish between the 
different reconstructions. 

To investigate which reconstructed velocity fields are 
consistent with regularly rotating gas discs, we tried to fit 
each of them with tilted ring models. In these fits, gas is 
assumed to move on circular orbits around the centre of the 
galaxy in a series of concentric rings. The position angle 



10 Ole M oiler & Edo Noordermeer 



o 

0) 

m 
o 
u 




X (arcsec) 



Figure 10. Residuals between the reconstructed velocity fields and tilted ring model fits. The panels are for the same reconstructions 
as in the previous figure. The lines show the isocontours of the reconstructed source image at 1, 3 and 5 scale lengths. 



and inclination of the galaxy, as well as the rotation ve- 
locity of each ring is fitted to obtain the best match with 
our 'observed' reconstructed velocity field. The differences 
between the lens reconstructions of the source kinematics 
and our best tilted ring fits are shown in Fig.Q^lfor all re- 
constructions. For the majority of cases, the differences in 
the central regions of the source are small, but larger in the 
outer, high- velocity regions of the source. However, several 
reconstructions (e.g. 2, 5 and 8) show significant distortions 
in the inner regions. Such distortions are not observed in 
real galaxies and must therefore be due to an imperfect lens 
model. 



5.3 Comparison of isophotal and kinematic fits 

Even though the source light profile is acceptable for all re- 
constructions, close inspection of the velocity fields shows 
that some of the lens models are insufficiently accurate to 
allow reconstruction of the source kinematics. To investigate 
this further, we also performed an isophotal analysis of the 
reconstructed images and compared the morphological ori- 
entation of the reconstructed sources with the kinematical 
orientation as derived from the tilted ring fits. The derived 
values for the position angle and inclination from both anal- 
yses are listed in Table 2. Comparing the position angles of 
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Table 2. Parameters of the source reconstructions for the 9 
best-fitting lens models. The values for the position angle from 
isophotal and tilted ring fits, 4>iso an d </>kin are listed together with 
the corresponding inclination angles, #i so and #kin- Also shown are 
the results for the substructured models CI, C2 and C3, discussed 
in § 5.5. 



the isophotal analysis to those obtained from the tilted ring 
fits shows a large discrepancy of 10 or more degrees for sev- 
eral reconstructions (2, 5, 6). The inclination angles from 
the isophotal and tilted ring fits agree to within 10 degrees 
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Figure 11. Rotation curves from tilted ring fits. The thin black 
lines give the rotation curves derived from the velocity fields from 
reconstructions 3, 4, 7 and 9 from Table 2. The bold line with 
data points shows the average of the 4 individual curves; the 
errorbars are estimated as 1/2 of the rms scatter of the points 
from the individual curves. The solid, thin green line shows the 
input rotation curve from ea. llll The red line shows the result for 
the average of the three best fits to the substructured source (cf. 

in all cases. The inclination angles from isophotal fits are 
within 10 degrees of the input value of 50 degrees with the 
only exception of reconstruction 2 which has a best fit incli- 
nation angle of 35 degrees. Only for models 3, 4, 7 and 9 is 
|#iso — #kin| < 5 degrees. Inspecting the residuals between the 
velocity fields as fitted with the tilted ring models and the 
reconstructed velocity fields, shown in Fig. 1101 also helps to 
discriminate different models. The strongest residuals in the 
central parts - that is inside the contours in each panel of 
Fig. llQI - are associated with the reconstructions for models 
1, 2, 5 and 8. 

On the basis of the information from Table 2 and 
Fig-EH one would draw the conclusion that reconstruc- 
tions 3, 4, 7 and 9 are the most accurate. From these, one 
would constrain the position angle and inclination angle 
of the source to be |(0kin + </>iso) = 0.2 ± 0.5 degrees and 
f (#kin + #iso) = 47 =b 3 degrees - consistent with the input 
values. In Fig.[H] we show the rotation curves derived from 
the tilted ring fits for these reconstructions, together with 
the average of all 4 and the input rotation curve from ea.lTTl 
Within the expected observational uncertainties these re- 
constructed velocity fields allow an accurate recovery of the 
input rotation curve. 

5.4 Effect of cluster mass profile and ellipticity 

Thus far we have assumed a cored pseudo-isothermal model 
for the cluster halo mass distribution in our model. The total 
mass profile in our model cluster follows an NFW form only 
due to added substructure in the form of galaxies. A pseudo- 
isothermal mass profile has been shown to provide good fits 
to lensing cluster halos in previous studies (|Sand et aI 2005: 
Kne ib et al.l lT996 > ) . Numerical simulations, however, predict 
that cluster halo mass profiles follow a NFW form more 
closely. In addition, we have assumed a spherical mass dis- 
tribution for the cluster, whereas most real and simulated 
clusters have elliptical mass profiles. There is a degeneracy 
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Figure 12. The reconstructed source kinematics of the best-fit 
elliptical NFW cluster lens model. The original simulated data 
(shown in Fig.0 is fitted with an elliptical NFW halo model 
with three cluster galaxies. The figure shows the reconstructed 
integrated light profile (contours) and velocity field (colour) of 
the source. The best-fit NFW lens model has an ellipticity of 
e = ^\-b 2 la 2 = 0.15, where a and b are the major and minor 
axes, respectively. Its virial mass M v [ r = 1.925 x 10 14 M(T) and its 
concentration parameter is c = 5.6 (see Nav arro et alJ Jl997f) for 
a definition of these parameters). 

here between the ellipticity of the halo mass profile and the 
presence of massive substructure, as both induce an asym- 
metric lensing potential. The effect of triaxiality on lensing 
statistics has also been found to be degenerate with core size 
bv lQguri et, all (2005V We have tested how strongly elliptic- 
ity in the lens model affects the properties of the arcs and 
the reconstructed source. Fig. 1121 shows the reconstructed 
source properties for the best fitting model of an elliptical 
NFW cluster halo with three additional cluster galaxies. A 
relatively good fit to the source light profile is retained. How- 
ever, the model fails in reproducing the small scale details 
of the velocity structure of the source; the velocity field in 
the central region is bend into an 'S '-shaped structure. In 
principle, high quality kinematic data could break the de- 
generacy between halo ellipticity and massive substructure. 
In practice, however, it is doubtful whether current instru- 
ments could provide data that would be accurate enough to 
detect these small scale differences. 



5.5 Fits to substructured sources 

Above we discussed how the appearance of the velocity field 
of arcs would change if sources are substructured instead of 
smooth. Since such substructures would produce additional 
uncertainties in any source reconstructions, we also perform 
fits to a substructured source model using our method. The 
velocity fields and residuals to tilted ring fits for the three 
best fitting reconstructions C1-C3 are shown in the bottom 
row of Figs.l^lfeHUl respectively. Models CI and C2 show a 
low residual to the tilted ring fit, whereas model C3 shows 
some velocity distortions in the central region and appears 
slightly more compact overall. Averaging the resulting best- 
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fit rotation curves for these three reconstruction gives an 
averaged inclination angle of 42 zb 4deg and an average po- 
sition angle of 4 ± 1 deg. These values are within 2a of the 
input. The averaged rotation curve is shown as the red line 

in Fjg.rrn 

5.6 The lens mass reconstructions 

We investigated what differences in the lens mass recon- 
structions lead to the observed differences in source recon- 
structions. Inspecting the reconstructed mass maps and cal- 
culating the total mass inside the Einstein radius for all 
the reconstructions, we found that there are only very small 
differences between the different reconstructions; the total 
mass within an Einstein radius is always within 1-3% of the 
input value. However, there are larger differences in how the 
mass is distributed within the Einstein ring; the reconstruc- 
tions differ in their amount and direction of asymmetry in 
the central regions. In addition, the surface mass density 
in regions away from the main arc vary by ~ 30 — 50% 
between the reconstructions. This large variation is to be 
expected as the main constraints on the lensing mass distri- 
bution comes from the main arc itself and is, strictly, only 
a local constraint on the mass distribution. It is only due to 
the parameterised form of the input model that the arc con- 
strains other regions in the cluster at all. The differences in 
the reconstructions of the velocity fields mainly come from 
the change in the distribution of mass in the central part 
of the cluster and from small changes in the mass of the 
galaxies very close to the main arc. 



6 OBSERVATIONAL POSSIBILITIES 

In the previous sections, we have shown how gravitational 
lensing affects the velocity fields of high-redshift background 
galaxies, and how the additional information contained in 
the observed velocity fields can be used to constrain models 
of the lensing cluster and source galaxy. In this section, we 
will discuss the possibilities of observing the velocity fields 
of strongly lensed galaxies. 

It is clear from the previous sections that one needs to 
measure velocities over the full 2D extent of an arc, in or- 
der to use the kinematic information to constrain the lens 
model and to determine the rotation curve of the source 
galaxy. Simple long slit spectra along lensed arcs lack infor- 
mation about the orientation of the source galaxy, making 
it impossible to interpret the kinematic properties of the 
arc unambiguously. Additionally, small slit offsets and finite 
slit widths have a strong effect on the observed profiles, as 
can be appreciated from the difference between for example 
panels (a), (c) and (e) in Fig.|H| 

To observe the velocity fields of giant arcs, integral field 
spectroscopy at high spatial resolution and high sensitiv- 
ity is required. Our simulated lensed velocity maps in Fig-El 
have a width of about 5 arcseconds, but one should bear in 
mind that in producing these figures, we have not applied 
any flux-cut. In reality, velocities can only be measured in 
the brighter regions of the arcs, with a typical width of 1 
arcsecond. Spatial resolution of at least 0.2-0.3 arcseconds 
is required to resolve these regions. The demands on spec- 
tral resolution are less stringent. Typical spiral galaxies have 



rotation velocities in the range 100-300 km/s, so a velocity 
resolution of order 50-100 km/s is sufficient to measure ra- 
dial velocities of the gas to a fraction of the expected rotation 
velocities. 

To carry out the required observations, a number of op- 
tions are available. Currently, the best opportunity is offered 
at optical or near-infrared wavelengths, where, depending on 
the source redshift, several strong emission lines of ionised 
gas (e.g. Ha, On, Om, etc.) are available. Sub-arcsecond see- 
ing is now routinely achieved with adaptive optics systems 
at a number of ground-based telescopes, and the number of 
integral field spectrographs that make use of the high res- 
olution offered by these systems is rapidly increasing (e.g. 
GMOS on Gemini and Sinfoni on the VLT). The biggest 
obstacle currently seems to be that, in order for adaptive 
optics systems to deliver the sub-arcsecond images, a bright 
(mv < 14) guidestar in the immediate neighbourhood of the 
object is required. Since most lensed arcs do not lie close 
enough to such a bright star, one has to await the devel- 
opment of artificial laser guide-star systems to observe the 
most interesting arcs. However, technology seems to be im- 
proving rapidly, and several observatories expect a working 
system within a few years from now. 

All giants arcs observed hitherto are intrinsically faint, 
so large telescopes are required to obtain useful spectra. 
Long- slit sp ectra of a straigh t arc at z = 0.91 have been 
obtained bv lPello et alJ (|l99lh . using a 2arcsec wide slit at 
the 4.2m WHT. In 15 hours of integration time, they ob- 
tained a high signal-to-noise ratio spectrum which enabled 
them to extract velocities along the full length (> lOarcsec) 
of the arc. Several other groups have recently measured spa- 
tially resolved velocities in unlensed galaxies out to red- 
shifts of z rsj 2, using 8-1 0m class telescopes like Keck or 
VLT and slitwidths of 0.5-1.0 arcsec (IVogt et alll99fl [l997 ; 
iBohm et al1l2004t lErb et al.ll2004l) . These results imply that 
at 8-1 0m class telescope like the VLT, Keck or Gemini, sub- 
arcsecond resolution observations should be feasible in 1-2 
nights of integration time. 

Other options lie further in the future. ALMA is cur- 
rently being constructed, and will offer the required resolu- 
tion and sensitivity to observe the kinematics of molecular 
gas at the redshift s of the arcs we study here. Even further 
ahead, giant radiotelescopes like SKA will be able to observe 
the Hi emission line of neutral hydrogen. This would offer 
the fascinating possibility of measuring the kinematics of 
lensed galaxies well outside their stellar discs, probing into 
the dark matter dominated regions of these young galax- 
ies. Finally, several studies are currently underway to design 
the next generation optical telescopes, with diameters of 25 
meter and larger. With the light gathering powers of such 
extremely large telescopes, it will be feasible to detect emis- 
sion lines out to large galactocentric distances in lensed arcs 
within very short exposure times, thus enabling systematic 
studies of the kinematics of these high-redshift galaxies. 



7 DISCUSSION AND CONCLUSIONS 

Determining the properties of high redshift galaxies remains 
one of the main goals of current research. In this paper, we 
presented a first theoretical investigation on how the effect 
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of gravitational lensing can be exploited to determine the 
kinematic properties of high redshift galaxies. 

Using a parametric cluster model we simulated the ve- 
locity structure of strongly lensed background galaxies. The 
combination of ray-tracing techniques with parametric clus- 
ter and source models proved to be a very efficient and ac- 
curate method for this study. 

In general, we found that the two dimensional kine- 
matic profile along strongly lensed arcs is very complex. 
Differential magnification leads to very distorted position- 
velocity diagrams and strong asymmetries in the velocity 
fields. Here, it is important to note that the individual clus- 
ter galaxies close to the arc contribute strongly to this effect, 
as we demonstrated in sections 4 and 5. Using a relatively 
simple-minded technique, we showed that reconstructions 
of the 2D kinematic source properties of lensed arcs are in 
principle possible. Since the velocity structure is sensitive to 
small variations in the lensing potential, kinematic informa- 
tion along the arc provides additional tight constraints on 
the mass distribution in the proximity of lensed arcs. Ob- 
servationally, the use of an integral field spectrograph at an 
8-10m class telescope with sub-arcsecond seeing will allow 
accurate source reconstructions and measurements of the 
rotation curve of strongly lensed arcs. We predict that the 
inclination and position angles of sources that are dynam- 
ically stable rotators can be determined to an accuracy of 
~ 10%. Rotation curves can thereby be determined with ac- 
curacies of better than 15% out to 2-3 disc scale lengths for 
galaxies at redshifts above z ~ 1.5 in this way. 

Our general approach was motivated by our aim to 
provide a general discussion of the kinematic properties of 
lensed arcs and point out the uses and possibilities of kine- 
matic data of such systems. The parametric cluster model we 
used was simple, but it reproduced the observed appearance 
of strongly lensed arcs well, when the galaxy mass-fraction 
inside the cluster Einstein radius was high. We note here 
that it should be possible to constrain the mass fraction of 
galaxies in clusters by making statistical predictions about 
the shape of arcs, for example from N-body simulations, 
and comparing them with the observed structure of arcs. A 
thorough study of this would have to take into account the 
effect of cluster merging and the mass function of cluster 
sub-haloes. Our predictions for the general appearance of 
velocity fields of arcs are independent of the specific form of 
the cluster potential, as long as the observed properties of 
arcs are reproduced. This is because any other description 
of the cluster potential must also reproduce the observed 
properties of lensed arcs. In particular, the local differential 
magnification that produces the distortions in the velocity 
fields arises from small scale mass structures close to the 
arcs, which are also responsible for the broken structure of 
observed arcs. 

Our method here has made use of a smooth source 
model. Even though we demonstrated that structure in the 
source does not change our results and does not affect the 
appearance of arcs and their velocity fields significantly, this 
is strictly only true as long as the source itself has stable 
rotation. For high redshift sources that will probably not 
always be true, since mergers, outflows etc. are much more 
common at higher redshifts. However, as we pointed out, 
one of the advantages of the multiple arc systems formed 
by lensing is that intrinsic properties of the source can be 



disentangled from lensing induced distortions. Lensing dis- 
tortions will be different for each arc, since the local cluster 
mass structure is important, whereas intrinsic source prop- 
erties are the same for all arcs from a single source. In fact, 
this can be exploited to the extent that the source and lens 
can be reconstructed from arcs without (almos t ) an y as- 
sumptions about the source itself. IWarren fc pTel fcOpA de- 
scribe a method that exploits this and can reconstruct the 
lens and source in the presence of noise and finite seeing. 
Such a non-parametric metho d was previously d escribed by 
Wall ington et al.1 (|l996L I1994T) and extended by iKoopmansI 
( 2005). It can be applied equally well to reconstruct the kine- 
matic source properties, independent of any assumptions 
about the source - except that the source be of a physically 
plausible size. 

In this paper we concentrated entirely on arcs produced 
by clusters. This was motivated by the fact that cluster arcs 
are larger and easier to distinguish from light emission origi- 
nating from the lens plane. Galaxy lenses produce consider- 
ably smaller arcs that are superimposed on the lens galaxy 
itself. However, with IFU's of high spatial resolution it may 
be possible to determine the kinematic profile of arcs lensed 
by galaxies as well. Since the relative scale of source to lens 
is about unity for galaxy lenses, in contrast to cluster lenses 
where it is much smaller, the appearance of the arcs pro- 
duced by galaxy lenses is much smoother than for cluster 
lenses. This can be explained by noticing that a version of, 
for example, the top left panel in Fig. 0] that is scaled down 
by a factor of 10 or more would be covered almost com- 
pletely by the source. This means that small scale structure 
in the lensing potential would have almost no effect on the 
overall appearance of the lensed arcs. However, if kinematic 
data is available, the situation changes. Differential magnifi- 
cation becomes important for each individual channel since a 
given velocity channel only probes a very small region on the 
source plane. In other words, kinematic data of arcs behind 
galaxy lenses can provide strong constraints on the amount 
of mass in small scale structures of galaxy haloes, and may 
be used as a direct probe of the mass function at the low- 
mass end - possibly down to ~ 1O 7 M . We will address this 
in more detail in a future publication. 

In summary, it is clear that obtaining two dimensional 
kinematic profiles of strongly lensed arcs will provide very 
useful information of both source and lens. Using IFU's in 
the very near future to study these systems will provide 
a unique way to determine the rotation curve of strongly 
magnified galaxies at redshifts z ~ 1.5 or higher and measure 
their mass-to-light ratio out to several disc scale-lengths. 
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